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The specification of polishing requirements for the optics in coronagraphs 
dedicated to exo-planet detection requires careful and accurate optical mod- 
elling. Numerical representations of the propagation of aberrations through 
the system as well as simulations of the broadband wavefront compensation 
system using multiple DMs are critical when one devises an error budget for 
such a class of instruments. In this communication we introduce an analytical 
tool that serves this purpose for Phase Induced Amphtude Apodisation 
(PIAA) coronagraphs. We first start by deriving the analytical form of the 
propagation of a harmonic ripple through a PIAA unit. Using this result 
we derive the chromaticity of the field at any plane in the optical train of a 
telescope equipped with such a coronagraph. Finally we study the chromatic 
response of a sequential DM wavefront actuator correcting such a corrugated 
field and thus quantify the requirements on the manufacturing of PIAA 
mirrors. 

(c) 2011 Optical Society of America 
OCIS codes: 220.2560, 070.6760, 120.6085 
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1. Introduction 

Analytical propagation of wavefront errors through coronagraphs is the basis of all the recent 
error budget estimates for both ground and space based exo-planets imaging instruments. 
The literature has been particularly active in this area and our understanding of the sensi- 
tivity to aberrations of several coronagraphic solutions has considerably expanded over the 
past few years. For band limited coronagraphs, Sivaramakrishnan et al. [1] first introduced 
an analytical propagator for low-order modes based on a Zernike decomposition. In parallel, 
Shaklan & Green [2] carried out the same analysis for mid-spatial frequencies, based on 
a Fourier decomposition, comparing the sensitivities of band-limited coronagraphs. Perrin 
et al. [3] studied the impact of second order terms on the point spread function, a result 
revisited and formalized by Give'On et al. [4], who introduced the concept of frequency 
folding. Sivaramakrishnan et al. [5] then focused on the propagation of mid-spatial frequen- 
cies through Apodised Pupil Lyot Coronagraphs. The case of out of pupil optics and the 
resultant chromatic phase to amplitude mixing was first tackled by Shaklan et al. [6] for 
band-limited and shaped pupil coronagraphs and we recently provided a formal analytical 
approach that addresses this issue in Pueyo & Kasdin [7]. 

Computing the sensitivity to aberrations proves to be more challenging task in the case 
of Phase Induced Amphtude Apodized (PIAA) coronagraph and for the Optical Vectorial 
Vortex Coronagraph (OVVC). PIAA coronagraphic technology, first introduced by Guyon 
[8] , is a promising solution since it makes most of the photons collected by the primary mirror 
of the telescope available for planet detection and characterization. This technique is based 
on two aspherical mirrors that redistribute the light in the pupil plane of a telescope so that it 
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follows a given amplitude profile that leads to a Point Spread Function (PSF) having contrast 
levels close to 10~^°. Because all the light collected is remapped using these mirrors, such 
a coronagraph has virtually no throughput loss. As a consequence, the angular resolution 
is undiminished and is close to the diffraction limit of IX/ D, a feature comparable to the 
performances of a phase mask coronagraph such at the OVVC [9] , without any transmissive 
optics. The OVVC fully extinguishes on-axis starlight by introducing an azimuthal phase 
ramp at the focus of the coronagraph. This phase profile is obtained by manipulating the 
transverse polarization state of the hght with space- variant birefringent elements [9] . The 
difficulty in modeling OVVC resides in a proper treatment of the polarization effects and 
devising high fidelity models for manufacturing defects at the singularity located at the 
center of the focal plane mask. These issues are discussed and addressed in Mawet et. al 
[10]. 

The difficulty in PIAA modeling PIAA coronagraphs resides in devising accurate analyt- 
ical models for the propagation of low-order and mid-spatial frequency aberrations. Since 
the optical surfaces of such a system are highly aspherical, classical tools based on the 
Fresnel approximation cannot be used. The first numerical diffractive study of wavefront 
propagation through such a coronagraph was carried out by Belikov et al. [11], based on an 
expansion in Zernike polynomials. They showed that the high sensitivity to off-axis sources 
of pupil mappers was the cause of a higher sensitivity to low order aberrations. Herein we 
expand upon the results of Shaklan et al. [12] and derive a full treatment of the propaga- 
tion of mid-spatial frequency harmonic aberrations through a two mirror pupil remapping 
system. Our main result is an analytical propagator for harmonic wavefront errors entering 
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a PIAA coronagraph. The derivation of this analytical tool is presented in Section 3. The 
remainder of the paper illustrates how this propagator can be used to address the prob- 
lem of broadband wavefront control for pupil mapping coronagraphs. Section 4 follows the 
presentation of Pueyo & Kasdin [7] to derive an analytical expansion of the wavefront of a 
telescope equipped with such a coronagraph. Finally, in Section 5 we apply this analytical 
method to predict the broadband performances of wavefront sensing and control systems 
applied to a pupil mapping coronagraph. 

2. Background 

As shown in Fig. 1, a pupil remapping unit is composed of two aspherical mirrors that 
remap the light of an incoming pupil according to a prescribed apodisation profile Traub 
& Vanderbei [13]. As presented in Guyon [8] a PIAA coronagraph is composed of two 
remapping units separated by a focal plane mask. The first set of aspherical mirror remaps 
the telescope pupil so that the starhght is concentrated in the core of a very high contrast 
Point Spread Function. The focal plane stop blocks the core of this PSF and hence removes 
the bulk of the starlight while preserving most of the photons form a potential companion. 
The purpose of the second remapping unit is to invert the remapping in order to restore 
the imaging properties of the whole apparatus. The purpose of this paper is to investigate 
the physics and propagation properties of a single two-mirrors remapping unit and thus 
solely focus on the "forward" combination shown on Fig. 1. In a future communication we 
plan to use the findings reported here to study the performances of a full four mirror PIAA 
coronagraph. In this section we review the design equations of a circularly symmetric pupil 
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remapper for which the designed apodization is independent of azimuth, 

A{x,y) = A{f). (1) 

Here (x, y) and (f, 9) are the location of the rays at M2 in Cartesian and polar coordinates 
respectively. Similarly, (x, y) and (r, 6) are the location of rays at Ml in Cartesian and 
polar coordinates. In the most general shown in the bottom right panel of Fig. 1, 

the diffraction limited field at a location at M2, is given by the sum of diffracted 

wavelets at each point at the surface of Ml. In the particular case of ray optics there 
is no summation involved and there is only a one to one mapping between the field at 
(x, y) and the field at the incident point on Ml (xo(5, y), yo(5) y))- For ^ given the 
coordinates {xq[x, y),yo{x, y)) are derived from the mirror surfaces using Fermat's principle. 
(ro(5;, y), 6*0(5, y)) are the coordinates of the same point in a polar system. Without loss 
of generality we choose, for the remainder of this paper, to focus on pupil to pupil on-axis 
PIAA systems that are shown in the bottom two panels of Fig. 1. Moreover if we choose 
6*0 = 6*, then, following the presentation of Traub & Vanderbei [13], the relationship between 
the location of the incident and outgoing rays can be written, 

^ = ^A(ff. (2) 
or ro 

where ro is the radial location of the incident ray at Ml. The design of the mirror shapes is 
then driven by the following coupled partial differential equations, 

dh r(ro) — ro 



(3) 
(4) 
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dr ""^ Z 
dh . f — ro(r 

" z~ 



where h{r) is the height of Ml, /i(f) is the height of M2, and Z is the distance between the two 
mirrors. Eq. 3 and Eq. 4 were derived using a ray optics approximation. Unfortunately, in 
order to predict the contrast of an actual PIAA, it is necessary to carry out a full diffraction 
analysis. 

In particular, Vanderbei [14] showed that edge propagation effects limit the contrast to 
at most 10~^ for a two mirrors PIAA coronagraph alone. Pluzhnik et al. [15] subsequently 
showed that pre and post-apodisers can mitigate these effects and a 10~^° contrast can be 
recovered with little loss in throughput and angular resolution. Recently, we developed a 
new, purely analytical approach to the diffraction problem [16]. This new method can be 
applied to quantify the impact of diffraction of the edge of the remapping optics on contrast; 
using it we reproduced the results of Vanderbei [14] which established a set of pre and post- 
apodisers that allow 10~^° with an aberration free PIAA Pueyo et.al [17]. Here we do not 
delve into such an analysis, which will be the object of a future detailed communication, 
and choose to focus on the analytical modeling of wavefront propagation. We start from the 
main theoretical result of [16]. Namely, for any geometry, the diffraction limited field at M2 
is well approximated by the following quadratic integral, 

Eauti^, y) = 4^ I Ein{x, y)e'fz[i---o?+'^i---o)iy-yo)Hy-y.?]dx dy. (5) 

IaZ J Ml 

Where the ray optics remapping in cartesian coordinates is given by: 

xo(x,y) = ro(f)cos(e) (6) 
yo{x,y) = ro{f)sm{9) (7) 

The integral in Eq. 5 is a sum over all the contributions of each point at the surface of Ml, 
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while ray optics only relies on the field at {xo,yo)- Moreover, we show in App. A that for 
circularly symmetric systems the outgoing field is given by 

iXZ J Ml 

This result is the starting point for our analysis of the propagation of aberrations through 
PI A A systems. 

3. Diffraction analysis of PI A A 

In this section we are interested in finding an analytical expression for the outgoing field when 
Ein{x, y) is composed of harmonic ripples. While here we are considering the propagation of 
errors through a PIAA unit that is circularly symmetric we choose to describe the aberration 
in cartesian coordinates. By choosing a basis set that is not orthogonal over a circle we do 
not guarantee the unicity or the finite support of any given wavefront expansion. However 
this choice leads us to analytical insights about the wavefront mixing that occurs when a 
wavefront is propagated through a PIAA unit. Thus, here we follow the presentation of [6], 
and we study the propagation of harmonic ripples of the following form, 

^i^{Nr cos{4>-e)) ^g^ 

where N is the spatial frequency of the ripple and is its orientation. The propagation of 
such a complex disturbance through a circularly symmetric pupil mapping coronagraph is 
described by the following theorem. 

Theorem 1 Consider a pupil mapping system with prescribed apodisation A{f), whose in- 
cident illumination is harmonic: 

£;,„(r,^) = e^^(^'-'=°^('^-^» (10) 
8 



rA(r) 



\ (r cos{e-theta)-ro (f* sin{e-theta))'^ 



Then, assuming that the edge effects are mitigated according to the methodology presented 
in Pluzhnik et al. [15], using pre and post apodisers, the field distribution at M2 is: 

Where A{f) is the apodisation profile that is induced by the two remapping mirrors. 

Equation 11 combines the geometrical and diffractive propagation effects of pupil mapping 
systems for mid-spatial frequencies aberrations. This theorem is derived using Eq. 8 with 
Ein — e^^^'^^'^'^y\ our purpose being to decompose any incoming wavefront before the 
remapping mirrors into a sum of Fourier harmonics and propagate them analytically. 

Proof. We start with Eq. 8 and rewrite it into a form similar to the Presnel integral of 
harmonic aberrations. First we express the harmonic ripple in polar coordinates, 



-^(rcosid-0)-rof + ^{rsini9-9)f 



rdr dO 



°" ' iXZ Jo Jo 

iXZ J Ml 

where we have applied a series of coordinate rotations in order to write Eouti^i in a form 
convenient to complete the square in the integrand. We then extract all the terms that do 
not depend on {u, v) out of the integral and complete the squares. This yields 

E^t{f,9) = J_^i^^Nrocos(9-<i>)^-i^[cosH9-<p)r^+sin\9-<prf] ^^3^ 



iXZ 

X e 

J Ml 



^iu-r,-^t^.os{9-4>)? + i-^{v-^^-^sMe-4>)? 



du dv. 



The first exponential factor corresponds to the ray optics remapping of the ripple. The 
second exponential factor accounts for propagation induced phase to amplitude conversion 



and the integral captures the edge oscillation effects due to the propagation, in the sense 
presented by Pluzhnik et al. [15]. The propagation integral can be rewritten, 

iXZ J Ml J^x^ffi 

(14) 



where S-^ ^ q is an ellipse centered at (^^^-^^ cos(i9 — 0))^, sin(^ — 0))^) of semi major 

axis as — {jz fA{f)'^ )'^ semi minor axis hs — (^^)~^- When the phase oscillations at 
M2 have been mitigated using pre and post-apodisers, along the lines of Pluzhnik et al. [15], 
the diffractive properties are equivalent to a ray optics propagation, namely — >■ oo and 
ftg — >■ oo 

J_ f eH-^^''-'''^'^^Adu dv ~ A(r). (15) 
iXZ Jmi ^ ^ ^ ^ 

This finishes our proof. 

In the remainder of the paper we use Theorem 1 to compute the effects of propogated har- 
monic aberrations on the final image and to determine limits on the wavefront control system. 
Another application of Eq. 11 is presented in App. B where we evaluate the sensitivity of 
PIAA systems to off-axis sources. 

4. Propagation of harmonic aberrations 

In a previous paper [7] we illustrated the impact of out of pupil plane optic on the final 
aberrated field. These are responsible for mixing the amplitude/phase nature of aberrations 
and for changing their chromatic behavior. We showed that these effects could be mitigated 
using a two sequential DM wavefront controller. The same problem arises for PIAA corona- 
graphs, since it is intrinsically composed of two surfaces that are not only out of conjugacy 
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with respect to the telescope pupil, but that also are highly non-parabolic, thus introducing 
extra perturbations on the wavefront. In this section we show how to use the propagator 
derived in Eq. 11 to quantify the chromaticity of post-PIAA aberrated wavefronts. 

A. Phase to amplitude conversion 

We first consider the case of a phase error in the field right after Ml that is given by 
Ein{x,y) = e^^^^^^y) . Note that this phase error can either stem from optics before PIAA 
Ml or on the surface of Ml, 

h{x,y) = EAo&m,ne^^("^«) (16) 

m,n 

where 9m,n = Aq is the central wavelength and b_rn,-n = Kn,n dimensional 

Fourier coefficients. We assume for now that the coefficients are small enough that the 
field at Ml can be approximated by 

e^^'^~(l + i^/i). (18) 
A 

Note that under this linear approximation, a phase errors is equivalent to an imaginary 
disturbance of the pupil plane and an amplitude error to a real disturbance. Assuming a 
circularly symmetric PIAA, Eq. 8 yields. 



EUf,e) = ^J (l + 4Mx,.))e'*l^"~"-"-'»'"-^'"'""-"»%*d9. (19) 

iXZ J Ml A 

Therefore, using the Fourier expansion of 5h{x,y) in Eq. 17 and Thm. 1, 

Eoutif, e) = A{r){l + &m,ne'*^^^^^^^°"°'^' cos e^,„+ro sine sin0^,„) 

A 



m,n 



g-i^^^^^^^[(^)^COs2(e-er.,n) + (^)^sin^(e-e^,„)]^_ (20) 
11 



The general form for the propagator of phase errors in classical coronagraph behaves, as 
derived in Pueyo & Kasdin [7], as e . The propagator derived here does feature the 

same behavior but also captures the high curvature of the PIAA optics through the angular 
magnification factor (^^^^)^ cos^(6' — ^m,n) + (^)^sin^(^ — ^m,n)- This magnification is 
different in the radial and tangential directions. Therefore, in the angular spectrum factor, 
the contribution of the radial magnification is weighted by the relative orientation of the 
ripple, 9m,n, with respect to the line of observation, 9. This yields a term in cos^(^ — 9m,n)- 
The same consideration for the tangential direction yields a term in sin^(6' — 9m,nY- Seen 
from M2 it is as if the ripples at Ml were propagating along two orthogonal pupil mappers 
of different linear magnification laws. As a consequence, a phase ripple at Ml will not only 
appear at M2 as a condensed oscillatory pattern, as predicted by the laws of geometric 
optics, but will also see some of its energy transfered to amplitude, thus creating phase 
induced amplitude errors. This effect gets stronger with spatial frequency as illustrated in 
Figs. 2 and 3. In Fig. 3, computed at A = 700 nm, for a spatial frequency of 10 cycles per 
aperture, a pupil diameter of D = 3 cm and a mirror separation of 2; = 1 m, the effective 
propagation distance is equivalent to a quarter of a Talbot distance and the conversion 
is total: all the wavefront error becomes amplitude. This phase to amplitude conversion 
behaves as D"^ /\z, the Fresnel number of the PIAA unit. Thus this design parameter has a 
direct impact on the feasibility of broadband wavefront control. 
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B. Amplitude induced phase error 

The same approach can be carried out starting with amphtude errors in the field right after 
Ml, 

Einix, y) = ro(l + ^ a^,„e^^(--+"f)) (21) 

m,n 

with a-rn,-n — o^m^n dimcnsional coefficients and ro is the average transmissivity of 

the incident field. This leads to the derivation of amplitude induced phase errors that is 
given by 



Eout{f, e) = A(f)(l + a^,„e^^^"^'+"'(''°^°"^" cos 0^,„+ro sine sine^,„) 

1). (22) 



This amphtude to phase conversion is illustrated in Figs. 4 and 5, computed at A = 700 nm. 
Once again, for a spatial frequency of 10 cycles per aperture, a pupil diameter oi D — cm 
and a mirror separation of ^ = 1 m, we observe that for a point located at the center of 
M2, the effective propagation distance is equivalent to a quarter of a Talbot distance and 
the conversion is total: all the wavefront error becomes phase. These considerations raise a 
fundamental issue when one seeks to create a broadband null in the image plane of a PIAA 
coronagraph. Because this phase to amplitude mixing is chromatic, it alters the bandwidth 
of wavefront correctors. This is the question we address next by deriving a full expansion of 
the chromaticity of the wavefront after a PIAA coronagraph. 
C. Wavelength expansion of the propagated electrical field 

In the previous subsections we derived the field propagation for a harmonic aberration at 
a single wavelength. Here we seek an expansion for the propagated field over a band of 
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wavelengths. Our main result is stated in the following theorem. 

Theorem 2 We assume that the edge effects are mitigated in the PIAA unit via pre and 
post-apodisers and that the optics before PIAA Ml and after PIAA M2 are such that m n <^ 




, where D is the optic diameter and z is the distance between a given optics and the 



conjugate of the telescope pupil. Then, the field in any plane of the optical train of a telescope 
equipped with a PIAA coronagraph, after an arbitrary number of reflections on aberrated 
optics, and a propagation through the two mirrors remapping unit, can be expanded using 
the following X-Fourier expansion 

E{x,y) = Aix,y) (l + ^ J2^'^_^^^^^im^+ny)\ (23) 

where fz^_^ = {fm%)*- That is, the odd terms in the wavelength expansion are imaginary 
and the even terms are real. 

Note that this theorem is exactly the same as Theorem 2 in Pueyo & Kasdin [7] . It is as if for 
wavefront mixing purposes, a PIAA coronagraph behaves analytically exactly like a classical 
coronagraph; except that the chromatic phase/amplitude conversion due to propagation is 
magnified by a factor of Ai, the angular magnification of the PIAA unit. This effect will be 
quantified in the next section. 

Proof. In Pueyo & Kasdin [7] we established, using an induction argument and assuming 
that the optical surfaces are all parabolic or flat, that Eq. 23 was true in any plane of a 
classical coronagraph as long as m,n <^ . Here we are interrested in proving that 
the propagation through the non-parabolic optics of a PIAA coronagraph conserves this 
property: namely if the field E^n at Ml is such that Eq. 23 is true, then the field Eo^t after 
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M2 also satisfies this property. We write the field at Ml as: 

fMl,-k\k 

EUx,y) = 1 + E E^' 'u " 6^^^"^^+"^^ (24) 

m,m k 

As a result of theorem 1, the field at M2 is: 



\k 

m,m k 



£Ml,-kyk 

^ — J^i^f'^^ -\- y~^^fc -^m,n ^i^jV^,„(T-o cosg cos </)m.n+ro sing s\n(j)m.n) 



m,n k 



We write the propagator at a spatial frequency (m, n) as: 

If we expand this exponential in a Taylor series then Eq. 25 becomes: 



oo sMl-kyk-p 

E^tir, e) = A(f)[l + E E E i"'" '\k-p e'^^-"^'-"-"-^" -s^^.^+^osin^- sin0.,„)^^^^(^^ 

m,n k p=0 P-^ 

oo fMl,-k 

Eautif: e) = A{f)[l + E E ^''"^tIz^ E "''^ ^m,«(r, ^)Pe^^^'"'"(''° COS </>^.„+r-o sine sin^^.„)| 

(27) 

We then re- write as a new Fourier series the function of (f, 9), that is represented by the 
sum over (m, n) on the right of the factor. Since we have the freedom to arbitrarily 
index these new Fourier coefficients, for clarity we choose to call them f^i^J^'^ — fm'^^^^^~^^ 

f-Ml-k 

J m,n ^i^Nm.njrocosd cos 4'm.n+ro sin 9 sinipm.n)^, (^ 
m,n P^' 

— ^ y*^2,-(fc-p)^i^Ar^,^^,(f cosfii cos(/>^, „/+fsin6i sin</>^, (2g) 
m',n' 

The /jJV^''"^^ coefficients at M2 can be written as 

fMl-k 

fmi'^'^'"^ = Yl ^^^^^V'm,n(r, ^)f e*lT(™^o(i,2/)+n2/o(i,2/))g-i2j(m'S+n'y)^^ (29) 

m,n P- 
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where we chose to represent the integrand in cartesian coordinates for clarity. Finally we 
reduce the double sum over {k,p) to a single sum since their argument only depends on the 
difference k' — k — p. This yields 

/ rM2.-k\k' \ 

Eoutix: y) = A{x, 1 + E E ^ " e%K^+"'^) , (30) 

\ m' ,n' k' J 

which finishes the proof, since for optics after the PIAA unit we can use the results of 
Pueyo & Kasdin [7]. 



Note that Eq. 29 estabhshes an explicit relationship between the Fourier coefficients 
of electrical field distributions before and after a PIAA. Such a relationship is useful in 
wavefront control applications. As mentioned earlier, the only difference between classical 
coronagraphs and PIAA is that for a PIAA the errors are propagated through a modified 
angular spectrum. The result is that the phase to amplitude mixing depends on the location 
on M2 and is stronger/weaker than for a classical coronagraph by a factor of {^^-^^Y cos^{9 — 
0m,n)+(^)^ sin^(^— This location dependent wavefront mixing impacts the broadband 
performance of DM based wavefront controllers. It drives the size of the dark zone achievable 
using such controllers, which we discuss next. 
D. Largest correctable spatial frequency 

When the spatial frequency is small enough, the largest two terms in the expansion in 
Eq. 30 are ifm,m~^/\ phase errors, and f^^^, amplitude errors. However, when the spatial 
frequency of the aberration gets larger, the propagated terms become larger and higher 
order wavelength dependent terms grow. If two sequential DMs that follow the PIAA unit 
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are used to correct the phase errors, i /^^~^/A term, and the amphtude errors, fm,n^: ^ 
residual halo appears at high-spatial frequencies due to these higher order terms. Because 
of the N'^ dependence of the angular spectrum propagator this halo is a strongly increasing 
function of spatial frequency. In Pueyo & Kasdin [7] we defined the highest correctable 
spatial frequency, Adamite ; as the spatial frequency of a ripple for which the application of a 
two sequential DM correction no longer provides better broadband contrast than for the case 
without wavefront correction. Ripples above Numitc cannot be corrected over a broadband 
because the chromatic mixing of the wavefront is too large. For a classical coronagraph, 
designed using parabolic optics, this limiting spatial frequency can be written as follows. 

The same considerations are valid for a PI A A coronagraph. However, when the DMs are 
located after the remapping mirrors, then the average spatial frequency seen by the propaga- 
tor is increased by a factor of , as shown in Eq. 29 where the Fourier kernel is now written 
as a function of {xo{x,y),yo{x,y)). Thus, if the maximal correctable spatial frequency is 
expressed in terms of cycles per aperture before the PIAA, then NumitPiAA becomes 

^l^rmtPIAA=j;^^{^fl' (32) 

While the final wavefront expansion is similar for PIAA and classical coronagraphs, the 
higher order propagation terms are larger for PIAA, thus reducing the largest spatial fre- 
quency correctable under a broadband illumination. With DMs located before the PIAA 
unit the largest correctable spatial frequency would be driven by the surface errors at M2, 
back-propagated to the plane of Ml. We will study this configuration in a future commu- 
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nication. Here we emphasize that a two sequential DM wavefront controller that follows 
a PI A A unit will not be able to correct spatial frequencies above NumitPiAA, expressed in 
Eq. 32, over a broadband, due to a phase to amplitude mixing that is too strong. 

In the case studied here, a PIAA that does not include de-mapping mirrors after the 
focal plane stop [8], the image plane contribution of aberrations above NumUPiAA extends 
all the way to low spatial frequencies. This low spatial frequency leakage can be explained 
by using Fig. 2 to Fig. 5: at the edges of M2 the aberration appears mostly as a small spatial 
frequency ripple and consequently throws hght near the core of the PSF. While considerably 
damped by the apodisation profile at M2, such a leakage corresponds to a very chromatic 
wavefront that is highly uncorrectable under a broadband illumination using two DMs in 
series after the PIAA unit. In the case of classical coronagraphs, aberrations above Nnj^itC 
only have a small impact, due to the tail of the airy function, on the contrast in the Dark 
Hole of the coronagraph. In the case of a PIAA that does not include de-mapping mirrors 
after the focal plane stop, because of this low spatial frequency leak, aberrations above 
NiimitPiAA can potentially have an impact at low spatial frequencies and thus influence 
the performances of the wavefront control system. In the next section we present, as an 
illustration of the analytical propagator derived above, numerical simulations that quantify 
the effect of the chromatic wavefront mixing on the overall post-correction contrast in the 
following configuration: no de-mapping mirrors and a two sequential DM wavefront actuator 
that follows the PIAA unit. 
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5. Contrast predictions 

We use the analytical expansion in Eq. 23 to predict the best broadband contrast that can 
be achieved by a PIAA unit in the presence of fixed wavefront errors compensated by a pair 
of sequential DMs located after the remapping mirrors. Note that different, and potentially 
better, broadband performances can be obtained with DMs located before the remapping 
mirrors or a de-mapping unit after the focal plane mask. While these architectures can be 
studied using the approach presented here, their implementation requires novel wavefront 
control algorithms that are beyond the scope of this paper. Thus we decided here to focus on 
quantifying the limitations of the simplest solution possible. We will extend this study to all 
possible combinations of DM before / after the PIAA mirrors and with /without de- mapper 
in a future communication. 

Because the propagator is wavelength dependent, the wavelength expansion of the field at 
M2 exhibits an infinite number of term, whereas, as shown by [18], two sequential DMs can 
only correct for the A° and i l/\ terms. In this section we isolate one fourier component, at 
a given spatial frequency, and quantify how well a two sequential DMs wavefront controller 
can reject it under a broadband illumination. 

A. Methodology 

For a given phase error at Ml, and a bandwidth centered around Aq, we use a first order 
expansion of the wavefront 

= 1 + i^e^i?("^^+'^^\ (33) 
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We propagate the ripple according to Eq. 11. We repeat this process for each wavelength 
in the band considered and thus build a data cube of fields at M2 that is represented on 
the left panel of Fig. 6. Note that here with this first order expansion we are only studying 
the effect of PIAA propagation on small wavefront errors, leaving aside non-hnear effects 
due to wavefront excursions. We then compute the electrical field distribution at the final 
image plane for each wavelength. We are interested in an annulus between the IWA and 
the OWA in {\/ D)onSky, as represented in the right panel of Fig. 6, that is the region 
where we expect the wavefront control system to operate and create a Dark Hole. For 
the simulations shown here we chose IWA = 2 {X/ D)onSky and OWA = 7.5 {X/ D)onSky 
These angles correspond to angular separation on the sky and are related to the actual 
units at the science focal plane by (X/D) camera — •M.{X/D)onSky, where A4 is the angular 
magnification of the PIAA unit. Because the PSF of single spatial frequency ripples that 
have been propagated through a PIAA is extended, as shown on Fig. 11, a ripple such that 
y/m^ + > OWA will still leak in the Dark Hole. The calculation we are carrying out in 
this section quantifies how much this leak is correctable using two DMs in series after the 
PIAA unit. To do so we proceed as follows. Assume that the field at a given pixel of the 
image plane is rj) = E^'^{^, rj) + iE{"^{^, rj). Then we assume that a perfect dual 

DM wavefront controller affects this field in the following fashion: 

That is, we are assuming that the ideal wavefront controller cancels the electrical field at the 
central wavelength and features an achromatic leak for the real terms and a 1/ A leak for the 
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imaginary terms. Another possibility would be to assume a controller that features the same 
chromatic dependance but that minimizes the chromatic residual intensity over the entire 
spectral bandwidth. The outcomes of both approaches are similar in terms of broadband 
performances and thus here we focus on the one described by Eq. 34, which is illustrated in 
Fig. 7. In this example the optical design is such that D — 3 cm, 2; = 1 m, A = 700 nm and 
the spatial frequency N = 7 cycles per aperture. These values lead to a chromatic mixing 
of the wavefront that is so strong that even after an ideal wavefront controller, the contrast 
under a 10 percent bandwidth is still 10~^ at its worse. Since we have assumed a ripple of 
amplitude 1 we can conclude that, for the PIAA design used in the example of Fig. 7, phase 
errors of spatial frequency 7 cycles per aperture before Ml cannot create speckles that are 
larger than 10~^ in order to be corrected by a dual DMs wavefront controller located after 
the PIAA unit. Next we repeat this approach for a variety of spatial frequencies, optics size 
and separations and bandwidths. 
B. Results 

The first parameter studied here is the speckle extinction as a function of spatial frequency 
and Fresnel number. Because propagation effects scale with the Fresnel number, J' — 
we expect the brodband residual halo due to higher order terms to increase when the Fresnel 
number decreases. Fig. 8 illustrates this feature. It shows how the maximum of a composite 
PSF over several wavelengths, after the application of an ideal wavefront control, behaves 
as a function of spatial frequency and Fresnel number. For this figure we used a dark zone 
going from IWA = 2 {X/D)onSky to OWA = 7.5 {X/D)onSky This OWA correspond to 
the 32{X/D) Camera outcr limit due to the limited number of degree of freedom of a DM 

21 



with 64 actuators across the pupil, divided by the magnification, Ai, of the pupil mapping 
unit considered here. Fresnel number is a crucial parameter when designing PIAA units: 
as seen on Fig. 8, for T — 11250 a two sequential DMs wavefront controller manages to 
extinguish speckles by seven to eight orders of magnitude over a broad range of mid-spatial 
frequencies. Such a level of extinction coupled with reasonably small wavefront errors to 
start with, makes broadband wavefront control over a dark hole in the image plane with a 
PIAA coronagraph feasible. However for T — 140, this extinction is reduced to two or three 
orders of magnitude, which considerably hampers prospects for broadband wavefront control 
with such a design. The second parameter studied is the speckle extinction as a function 
of spatial frequency and bandwidth. Fig. 9 plots the worst contrast in the same dark zone 
as a function of the spatial frequency of the input ripple with each curve corresponding 
to bandwidths of AA/A — 0.1, 0.2, 0.3. In this example, for — 11250, the bandwidth 
does not have much of an effect on broadband contrast. However, as seen on Fig. 10, with 
T — 1250 the propagation effects of high spatial frequencies aberrations through the PIAA 
coronagraph become very large, and, as expected, the performance of the 2 DM controller 
becomes sensitive to bandwidth. 

This analysis provides a methodology for deriving prescriptions on the optics that precede 
a PIAA coronagraph, when the DMs are located after the pupil remapping unit. Our main 
conclusion is that for Fresnel number larger than 10"^, chromatic propagation effects do not 
have much overall impact on the effectiveness of a two DM wavefront controller. Indeed, as 
shown in Fig. 9, such an apparatus manages to cancel a given spatial frequency by 8 orders 
of magnitude over several bandwidths ranging from ten to thirty percent. In this regime 
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the broadband halo is more hkely to be dominated by non-hnear effects such as frequency 
folding, that scales as 1/A^, and is only partially correctable using two sequential DMs. This 
effect has been discussed in previous communications Give' On et al. [4] and they have not 
been taken into account in the present paper. Nevertheless, since propagation effects become 
larger when T get smaller, this study shows that in the regime of ~ 100, chromaticity of 
the residual halo is driven by the propagated wavefront errors from Ml to M2. For instance, 
for a PIAA designed with T — 140, a 10~^° contrast over a 20 percent bandwidth seems 
an unachievable goal with DMs located after the PIAA unit. Indeed, as shown on Fig. 10, 
the pre-wavefront control speckles due to ripples of 10 to 50 cycles per aperture before Ml 
would need to be small enough so that the raw contrast is lower than 10~^ in order to be 
correctable to the 10~^° level with a post PIAA dual DM wavefront controller. 

6. Closing remcirks and future work 

In this paper we derived an analytical propagator for aberrations through a PIAA coro- 
nagraph. This propagator, is based on a Fourier expansion, and captures the wavelength 
dependence of the field after the coronagraph. This aspect is of critical importance with re- 
spect to the design of ongoing and future experiments based on this coronagraph. It provides 
a technique for tolerancing and error budgeting coronagraphic optics when PIAA is used in 
conjunction with a broadband wavefront controller. In the case of two DMs located after 
the two remapping mirrors, we applied a contrast evaluation procedure which predicts that 
broadband wavefront control for this architecture is only possible for PIAA Fresnel numbers 
that are larger than 2000. In the near future we will explore algorithms that control two 
sequential DMs that are located before Ml, and study the case of a PIAA coronagraph 
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with de-mapping mirrors, such architectures potentially being more favorable and allowing 
smaller PIAA Fresnel numbers. 

Appendix A: Second order expansion of the PIAA integral with circular sym- 
metry 

We start from Eq. 5 and write the expression of the partial derivatives of the remapping, 
Eq. 2, in polar coordinates, 

dx 9f r . On 
— — — cos Uo H stn Uq 



dxo dro ro 

1^ = ^sinX + ^cos% (Al) 
oyo dro ro 

-K— = (t^ )cos6o sinOo 

oyo OTo ro 

where we have assumed that 9 = 9o- Eqs. Al corresponds to the partial derivatives of the 
inverse remapping and are derived using the chain rule. We are interested in changing the 
integration variables in Eq. 5 from cartesian to polar coordinates: 

{x — Xo)"^ — r^cos^O — 2rrocos9 cos9o + TqCos^^q 

{y ~ yof — r'^sin'^9 — 2rrosin9 sin9o + rlsin'^9o (A2) 
{x — Xo){y — Vo) — r^cos9 sin9 — rro{sin9 cos9o + cos9 sin9o) + TqCos^q sin9o 
dx dy — rdr d9 

where the integral is taken over a circle of radius R and centered at the center of this circle. 
We will call such a domain of integration C^^y After some algebraic manipulations we 
can find a simple expression for the radial terms in the exponential factor of the quadratic 
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expansion: 



: ^cos'(^-^o) + -sin2(^-^o) 
dro ro 

dv 

rro ■ 2^ — cos(^ — ^o) 
dro 



As a consequence the radial field distribution after M2 becomes: 

E^t{f) ^ —- e^n^^ -^0^ Mrdrd^ A4 

EUr) = / e'^i^^'"''"^''"^ncixdy (A5) 
zAZ „j 

The main insight of this expansion is the fact that the propagation between the two mirrors 
of such a PIAA system reduces to the integration over an equivalent elliptical aperture. 
The geometry of this ellipse varies with f, the location on M2. Qualitatively, based on an 
energy conservation argument, we already know that the local effective aperture size for 
propagation purposes is stretched by a factor of in the radial direction. The elliptical 

integral in Eq. 8 formally illustrates this intuitive result, which states that since the area 
of integration has to be A{fY, the effective local aperture size in the tangential direction, 
normal to the radial, has to shrink by a factor of ^ . 

Appendix B: Sensitivity to Off- Axis response 

The off-axis magnification of PIAA systems is the feature that makes such designs so ap- 
pealing to the exo-planet community since it is the source of their intrinsic extremely high 
angular resolution. It was first explained by Guyon using energy and area conservation argu- 
ments ([8]) and then formally derived by Traub and Vanderbei ([13]) using ray optics. Here 
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we evaluate this magnification using Eq. 11 in order to compute the sensitivity to off-axis 
sources. Consider an off- axis source illuminating a PI A A system with a wavefront tilted by 
an angle '')sky Assume that the resulting field distribution at M2 is Fourier transformed by 
an ideal lens. Call ^camera the angular location of the centroid of the point spread function 
of such a source traveling through the pupil mapping unit. Note that here we define the 
angular magnification using the centroid of the off-axis PSF and not its maximum. While 
this definition is less accurate it has the advantage of providing a number that does not 
vary with the angular separation of the off-axis source. We thus define M., the angular 
magnification of a PIAA unit, as: 

M = 1^^:11^ (Bl) 

ISky 

Then, given a two mirrors PIAA design, M. can be computed as: 



r 1^ AiffJ-^ cos2 e + — sin^ e fdfde (B2) 
ttR^ Jo Jo y ^^(^j '^o 

This result is a direct consequence of the ray optics remapping factor of Eq. 11. Without 
loss of generality we can assume that (p — 0. If we write the angular separation of the off-axis 
source with respect to the optical axis in units of A/D, then the terms corresponding to the 
geometric remapping in Eq. 11 are written as: 

For this calculation we leave out the equivalent angular spectrum factor since it will only 
change the phase of the companion. From Eq. Al we know that: 

^ = 4^ cos^ e + ^ sin^ e (B4) 

ox rA[r) tq 

^ = (4^^ - -)cos ^sin 9 (B5) 
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Thus, at a given point on the surface of M2, (f, 9), the local spatial frequency is given by 
the magnitude of the gradient of jsky xo{x, y) 



\ rA{r) ro 

In order to compute M formally, we are interested in finding the centroid of this extended 
PSF. More formally, the centroid of the planet PSF will be located at the barycenter of 
N{r, 0) weighted by ^(f)^, which gives: 



ICamera = {N{f, 0)) = ISky^ A{f)\ -^^COsW + -SlV?e fdfdd (B7) 

Trit^ Jo Jo V rA[r) tq 

When we use a PIAA unit that follows the 10~^° prolate profile for /1(f), this yields 
M. — 2.63. This is exactly the value found when measuring the angular magnification 
using simulations such as the one shown on Fig. 11, where we have computed the PSF of 
two off-axis sources, of respective angular separation 2 and 4 X/D, propagated through a 
PIAA unit. The centroids of these PSFs appear in the final image plane at 2 x 2.63 and 
4 X 2.63 X/D. Note that because the main contribution to this angular magnification comes 
from the center of M2, a legitimate approximation for this value is ~ ^4(0). A similar 
proof was presented by Guyon ([19]), assuming ray optics. This result is a fundamental 
property of PIAA systems and is the source of their high performance. Because of the full 
throughput, spatial frequencies are magnified and planets that are very close to their parent 
star can be observed. 
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Figure 1 

Setup of the problem and notations] Setup of the problem and notations. Top Left: Three 
dimensional representation of a pupil to pupil off- axis PIAA system. Top Right: Side view 
of the geometrical remapping in a pupil to pupil off-axis PIAA system. Bottom Left: Side 
view of the geometrical remapping in a pupil to pupil on-axis PIAA unit: This is the 
configuration that is studied in this communication. Bottom right: Side view of all 
the rays contributing to the diffractive field at a point of coordinates {x, y) at M2. The 
ray corresponding to the geometrical remapping, which has coordinates (xo(5, y), yo(5, y)) 
in the input plane, is highlighted. 
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Figure 2 

Phase induced amplitude errors through a PIAA system. Spatial frequency at Ml A'" = 6, 
with D = 3 cm and z = 1 
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Figure 3 

Phase induced amplitude errors through a PIAA system. Spatial frequency at Ml N — 10, 
with A = 700 nm D = 3 cm and z = 1 



34 



Figure 4 

Amplitude induced phase errors through a PI A A system. Spatial frequency at Ml A'" = 6, 
with A = 700 nm D = 3 cm and z = 1 
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Figure 5 

Amplitude induced phase errors through a PIAA system. Spatial frequency at Ml N — 10, 
with A = 700 nm D = 3 cm and z = 1 
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Figure 6 

Cartoon representation of the wavelength data cube at M2 -left panel-, and at the image 
plane -Right Panel. The transverse axis is a virtual cut across a wavelength cube. These 
cubes are obtained by stacking field distributions at M2 and the image plane for several 
wavelengths across the spectral bandwidth of interest. 
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Figure 7 

Top: Illustration of the fit through the wavelength cube for one pixel. Top Left: raw contrast 
at one pixel in the image plane. Top Right: contrast at one pixel in the image plane after 
a perfect two sequential DM wavefront correction. Bottom: Residual intensity in the dark 
zone after subtracting the two dominant terms of the wavelength expansion. Note that the 
PSF of a ripple propagated through PIAA is much more extended than in the case of a 
classical coronagraph. The chromaticity of the leakage close to the optical axis has been 
modified by the propagator that introduced a higher order wavelength dependence. This 
drives the best speckle extinction achievable over a broadband. N = 7, D = 3 cm, ^ = 1 m, 
AA/A = 0.1 
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Figure 8 

Maximum of the broadband halo in the dark hole created by two sequential DM wavefront 
controller as a function of Fresnel number and spatial frequency of the wavefront error. The 
top curve corresponds to the maximum of the non-corrected PSF in the dark hole: note that 
high spatial frequencies leak in the dark hole due to the spatial extent of the off-axis PIAA 
PSF. The other three curves show the maximum of the residual halo after correction for, 
from top to bottom, T = 140, 1250, 11250. 
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Figure 9 

Maximum of the broadband halo in the dark hole created by two sequential DM wavefront 
controller as a function of bandwidth and spatial frequency of the wavefront error. The top 
curve corresponds to the maximum of the non-corrected PSF in the dark hole: note that high 
spatial frequencies leak in the dark hole due to the spatial extent of the off-axis PIAA PSF. 
The other three curves show the maximum of the residual halo after correction for, from 
bottom to top, AA/A = 0.1, 0.2, 0.3. The Presnel number for the PIAA unit is = 11250. 
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Figure 10 

Maximum of the broadband halo in the dark hole created by two sequential DM wavefront 
controller as a function of bandwidth and spatial frequency of the wavefront error. The top 
curve corresponds to the maximum of the non-corrected PSF in the dark hole: note that 
high spatial frequencies leak in the dark hole due to the spatial extent of the off-axis PIAA 
PSF. The other three curves show the maximum of the residual halo after correction for, 
from bottom to top, AA/A = 0.1, 0.2, 0.3. The Presnel number for the PIAA unit is 1250 
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Figure 11 

PSF of two off axis sources tliat are separated by 2X/D and AX/ D from the star 
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Figures 
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Off-Axis PIAA Off-Axis PIAA, side view 




On-Axis PIAA, side view On-Axis PIAA, side view 

d=0 d=0 




Fig. 1. Setup of the problem and notations. Top Left: Three dimensional representation of a 
pupil to pupil off-axis PIAA system. Top Right: Side view of the geometrical remapping in 
a pupil to pupil off-axis PIAA system. Bottom Left: Side view of the geometrical remapping 
in a pupil to pupil on-axis PIAA unit: This is the configuration that is studied in this 
communication. Bottom right: Side view of all the rays contributing to the diffractive field 
at a point of coordinates (x, y) at M2. The ray corresponding to the geometrical remapping, 
which has coordinates {xo{x,y),yQ{x,y)) in the input plane, is highlighted. 
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x/D x/D 



2. Phase induced amplitude errors through a PI A A system. Spatial frequency at Ml 
6, with A = 700 nm D = 3 cm and z = 1 
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x/D x/D 



4. Amplitude induced phase errors through a PIAA system. Spatial frequency at Ml 
6, with A = 700 nm D = 3 cm and z = 1 
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5. Amplitude induced phase errors through a PIAA system. Spatial frequency at Ml 
10, with A = 700 nm D = ?> cm and z = 1 
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IWA = 2 A / D, OWA = 7.5 X I D 




Imaginary part of the field at I\/I2 Imaginary Part of the Image plane field 



Fig. 6. Cartoon representation of the wavelength data cube at M2 -left panel-, and at the 
image plane -Right Panel. The transverse axis is a virtual cut across a wavelength cube. 
These cubes are obtained by stacking field distributions at M2 and the image plane for 
several wavelengths across the spectral bandwidth of interest. 
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Fig. 7. Top: Illustration of the fit through the wavelength cube for one pixel. Top Left: raw 

contrast at one pixel in the image plane. Top Right: contrast at one pixel in the image plane 

after a perfect two sequential DM wavefront correction. Bottom: Residual intensity in the 

dark zone after subtracting the two dominant terms of the wavelength expansion. Note that 

the PSF of a ripple propagated through PIAA is much more extended than in the case of 

a classical coronagraph. The chromaticity of the leakage close to the optical axis has been 
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Fig. 8. Maximum of the broadband halo in the dark hole created by two sequential DM 
wavefront controller as a function of Fresnel number and spatial frequency of the wavefront 
error. The top curve corresponds to the maximum of the non-corrected PSF in the dark 
hole: note that high spatial frequencies leak in the dark hole due to the spatial extent of 
the off-axis PIAA PSF. The other three curves show the maximum of the residual halo after 
correction for, from top to bottom, T = 140, 1250, 11250. 
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Fig. 9. Maximum of the broadband halo in the dark hole created by two sequential DM 
wavefront controller as a function of bandwidth and spatial frequency of the wavefront error. 
The top curve corresponds to the maximum of the non-corrected PSF in the dark hole: note 
that high spatial frequencies leak in the dark hole due to the spatial extent of the off-axis 
PIAA PSF. The other three curves show the maximum of the residual halo after correction 
for, from bottom to top, AA/A = 0.1, 0.2, 0.3. The Fresnel number for the PIAA unit is 
= 11250 
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Fig. 10. Maximum of the broadband halo in the dark hole created by two sequential DM 
wavefront controller as a function of bandwidth and spatial frequency of the wavefront error. 
The top curve corresponds to the maximum of the non-corrected PSF in the dark hole: note 
that high spatial frequencies leak in the dark hole due to the spatial extent of the off-axis 
PIAA PSF. The other three curves show the maximum of the residual halo after correction 
for, from bottom to top, AA/A = 0.1, 0.2, 0.3. The Fresnel number for the PIAA unit is 
1250 
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Fig. 11. PSF of two off axis sources that are separated by 2\/ D and AX/D from the star. 
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